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Abstract 

The classical equation of motion of a Davydov model in a coherent state 
approximation is analyzed using the multiple scales method. An exponen- 
tially decaying long range interaction (Kac-Baker model) was included. In 
the first order, the dominant amplitude has to be a solution of the nonlin- 
ear Schrodinger equation (NLS). In the next order the second amplitude 
satisfies an inhomogeneous linearized NLS equation, the inhomogeneous 
term depending only on the dominant amplitude. In order to eliminate 
possible secular behaviour the dominant amplitude has to satisfy also the 
next equation in the NLS hierarchy (a complex modified KdV equation). 
When the second order derivative of the dispersion relation vanishes the 
scaling of the slow space variable has to be changed, and a generalized 
NLS equation with a third order derivative is found for the dominant am- 
plitude. As the coefficient of the third derivative is small a perturbational 
approach is used to discuss the equation. A complete solution is given when 
the dominant amplitude is the one-soliton solution of the NLS equation. 



1. Introduction 



Many quasi-one- dimensional molecular systems are very complicated struc- 
tures built from complexes of atoms - we call them "molecules" - connected 
by hydrogen bonds. As an example we mention the complicate structure of 
a-helix in protein. A very simple model consists in replacing the three cou- 
pled chains of the complex a-helix structure by only one chain of the form 
■ ■ -H — N — C = O — H — N — C = O ■ ■■ . Also only one of the intramolecular 
excitations - that corresponding to the amide I oscillation - is taken into account 
(further on they will be called vibrons). An acoustic phonon field describing the 
oscillations of the molecules along the chain is also introduced. As the amide I 
oscillation energy depends on the stretching of the adjacent hydrogen bond an 
anharmonic interaction between vibrons and phonons appears. 



It is now easy to write an Hamiltonian of Frohlich type for this very simplified 
model. For a system with a conserved number of vibrons the model was proposed 
by Davydov [|l[]-[f|] more than 25 years ago. Eliminating this constraint of fixed 
number of vibrons the model was extended by Takeno 0,0, @ at the beginning 
of the eighties. 

The Hamiltonian of Davydov's model is given by 

H = E EB£B n — ^ Jmn(B^B m + B+B n ) 

+ XX^Pn + \w{u n +l ~ U n f) + XE(««+1 - il n )B+B n . (1) 

Here B n ,B+ are the annihilation/creation operators for vibronic excitation of 
energy E in the cell n, J mn term takes into account the long range interaction 
between vibrons, u n ,Pn are the displacement operator and the corresponding con- 
jugate momentum of the n-th molecule, m is an effective mass of the molecule, w 
the elastic constant and % a coupling constant describing the nonlinear interac- 
tion between vibrons and phonons. As an example of the long range interaction 
between vibrons we shall consider an exponentially decreasing model (Kac-Baker 
model) 

Jmn = J ~z c ^ r = e ^ . (2) 

2r 

Since the pioneering paper of Sarker and Krumhansl || this model was intensively 
used by a series of authors to investigate the thermodynamic properties and the 
soliton characteristics in several nonlinear 1-D systems - [fT4}| . 

It is a general belief that a coherent state approximation is very suitable for 
describing extended localized states in such systems. Such an hypothesis was 
done by Davydov. His ansatz for the state vector is 

W) >=E^(*) B n^f-^E^(*)Pi-^i(*)%) l°>> ( 3 ) 

|0 > being the vacuum state both for vibrons and phonons. Here a n (t), (3 n (t), 7r n (t) 
are now c-numbers, time depending, which will be determined from a variational 
principle. Using the average value of H as Hamiltonian in the classical equations 
of motion one gets 

oo 

ifia n = Ea n - ^ Jp{a n + P + a n - P ) + x(Pn+i - Pn)a n 
P =i 

M(3 = w(p n+1 -2p n + p n ^) + X (\a n \ 2 - \a n ^\ 2 ). (4) 
In an adiabatic approximation (/?„ — > 0) from the second eq. (f|) we get 

Pn+l Pn \0"a\ , 

W 



2 



which introduced into the first gives 



iTid n — Ea n — ^ J p (a n + P + a n -p) ~ ~\ a n\ 2a n- (5) 



a typical self-trapping equation. 



2. Multiple scale approach 

The linearized problem admits plane wave solutions with the dispersion rela- 
tion 

huj(k) = E- J— -Y e'^cos klp = E- J— - ( smh 7 _A . ( 6 ) 

2r ~y 2r \ cos /17 — cos kl J 

A plane wave solution with constant amplitude exists even for the nonlinear 
equation (5), but with an amplitude depending dispersion relation (Stokes waves) 

hu(k) = e- j^-Y e ~ 1Pcoskl P- ^-M 2 - ( 7 ) 

It is well known that these solutions are unstable at small modulation of the 
amplitude (Benjamin Feir instability) [[HJ. To see this in our case we write 

a^a(l-A,(t))e« 

where A n {t) are satisfying the following system of linear equations 

ikl -ikl X 2 2 

Looking for plane wave solutions 

J^^ = ce i{vln-€lt) _|_ ^ e -i(uln-nt) 

the instability, Q 2 < 0, appears when 

which is always satisfied at small wave numbers v. 

In order to describe this process of amplitude modulation, the multiple scales 
approach (reductive perturbative method) ]TB[-||2D[ is used. The solution is writ- 
ten as an asymptotic expansion in a small parameter e. 

an = e mn^t)^A^t 2 M,-) (8) 
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where the amplitudes Aj depend on some slow variables £, t 2 , £ 3 , ... In order to 
see how these can be defined let us expand the dispersion relation around a point 
k (k = k + Ak) 



u(k) = cuq + uj\Ak + cu 2 Ak 2 + u) r , 



1 d n u(k). 



I k=ko 



n\ dk n 

and considers Ak = fie a . Then the phase of the plane wave can be written 

kin -ut= (k ln - u t) + fie a (ln - u^t) - /j 2 u 2 e 2a t - /j 3 u 3 e 3a t + ... (9) 

Two distinct situations can appear. In the first case W2 ^ (we shall call it "a 
normal situation") and the second when u 2 = 0. This last situation is referred 
as the "zero dispersion point" case (ZDP). In the normal situation, with a = 1, 
and the following definition of the slow variables 

£ = e{ln-v g t), T 2 = eH, r 3 = e\... (10) 

it is easily seen that a competition between nonlinearity, and will occur 

in order e 3 . Introducing these in (|j) we ask that the equation has to be satisfied 
in each order of e. In the first order in e the equation is linear and we obtain the 
dispersion relation (0). In the next order e 2 the equation is satisfied if the velocity 
v g is given by v g = u)\. In the order e 3 we found that the leading amplitude A\ 
has to satisfy the well known nonlinear Schrodinger equation (NLS) 

Recently [^TJ the analysis for Takeno's model was extended to the next order e 4 . 



Similarly, in the present nonhomogeneous linear equation satisfied by the 

next amplitude A 2 is obtained, namely 

.dA 2 d 2 A 2 , nl , i2 .> .dAi . d 3 A 1 

i^T + ^ + + 2|^i! 2 A) = + luiz-^-. (12) 

In the left hand side (lhs) we recognize the linearized NLS equation (1-NLS eq.), 
while the nonhomogeneity in the right hand side (rhs) contains only the domi- 
nant amplitude A±, solution of the NLS eq. In solving the equation ( Jl2| ) we are 
confronted with two distinct problems. Firstly we have to take care to eliminate 
any secular behaviours raised by the presence of the nonhomogeneity. They can 
appear from terms in the rhs which are members of the null space of the 1-NLS 
eq. As is well known the symmetries of the NLS eq. are solutions of the 1-NLS 
eq., so the dangerous terms in the rhs of (0) are to be found between them. 
Such a symmetry is easily identified in ([12]) namely 

^ = -(^ + 60^1 _1 (13) 
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and the possible secular behaviour is eliminated if the r 3 dependence of A\ is 
given by 

"^-+-3^ + 6^1 _j=0. (14) 

This is a complex modified KdV equation and is the second equation in the hier- 
archy associated to NLS eq. As all the equations in the hierarchy have the same 
spectral problem the T3 time dependence of A% can appear only in the initial 
positions and phases characterizing the solution A\ of the NLS eq. After this 
" renormalization" of the dominant amplitude A\ we remain with a nonhomoge- 
neous linear equation, free of secularities, for the second amplitude A 2 

dA 2 d 2 A 2 , . ,9 , s , t odA\ , , 

1 + U2 ~d^ + qi ~ 1 2 + 21 11 2) = - 6 ^3c|Ai| 2 ^. (15) 

A simple method to solve it when A\ is the one-soliton solution of NLS eq. will 
be given in the next section. More comments will be presented there. 

In the last years several papers have discussed the significance of the higher 
order approximations to a nonlinear dynamical problem, which in the lowest 



relevant order is described by a completely integrable equation |[22||-|[24||, and the 
role played by the next equations from the corresponding hierarchies in order to 
eliminate possible secular behaviours. Our results are in completely agreement 
with these discussions. 



3. ZDP case 

Let us discuss now the case of zero dispersion point. For models in which the 
carrier wave is of the form of a lattice plane wave such a point k exist always 
in the first Brillouin zone. It corresponds to the maximum group velocity of the 
localized excitation. For the dispersion relation (||) it is given by 

1 



cos k l = ~ \\/ cosh 2 7 + 8 — cosh 7 J . (16) 

The analysis presented in the previous paragraph has to be slightly modified. In 
the phase (§) the uo 2 term is missing. Then the next term has to be used to 
define the time variable r 2 and consequently we have to consider a — |. The 
slow variables are now given by 

£ = e*(ln-uit), r 2 = e 2 t. (17) 

In order e 3 we shall have a competition between the nonlinearity, the time deriva- 
tive and the third order spatial derivative t^ 1 - Also the second derivative 

^l 1 can contribute in the same order e 3 because it is multiplied by co 2 (k), which 
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expanded around ko brings an extra es dependence (u>2(k) = 3u 3 /j,es + ...). Col- 
lecting all these contributions of e 3 order we arrive at the following equation 
satisfied by the dominant amplitude A\ 

.dA x d 2 A, 2 . d 3 A, 

l ~^T + 3 ^3-^- + qW A, = iu 3 —. (18) 

Introducing dimensionless quantities 



X = ^=k^ T = u 3 k 3 r 2 , y = J-J?A 1 (19) 
V6^ V nu 3 k$w 



the equation fll8"|) becomes 



d<5/ ld 2 ^ . l2 d 3 ^ , , 

7 af + 2^ + l *'* = ^ (2 °) 

3 

where (3 = (6fj,)2. This is a NLS eq. perturbed with a third order derivative 
term. As (3 is small the eq. ( |2H| ) will be solved by a perturbation approach. We 
consider 

-> * + (21) 

where now \I/ is a solution of the unperturbed NLS eq. Also we introduce a slow 
time variable r = /3T. In first order in f3 we get 

+ -^^r + (2 * 2 <5tf + ^ 2 c^*) = -i— + i— — . (22) 
dT 2 dX 2 ' dr dX 3 K J 

This equation is exactly of the same form as (0). The possible secular behaviour 
is eliminated requiring that the r dependence of $ is given by the complex mKdV 
equation 

<9^ a 3 * , l9 <9^ 

-a7 + Wo +m ax-°- (23) 

We remain with a linear nonhomogeneous equation for namely 

^ + ^ + PI*| 2 » + * 2 «*) = -«l*| 2 i (24) 

which is free of secularities. 

Further we shall give a complete solution for the case when \1/ is the one- 
soliton solution of the NLS eq. If z — u + is the eigenvalue of the spectral 
problem the one soliton solution is given by 

— i<E> 

^ = 2v 6 - - , $ = 2wX + 2(w 2 -v 2 )T + $o, 6 = 2w(X-A a + 2mT). (25) 
cosh 9 
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As mentioned before, the r-dependence can appear only in the initial phase <3>o 
and the initial position X . Introducing ( p5|) in (|23"D it is easily find that 

<w 2 2\ dX . 2 2 
— — = Su(u —3v), — — = 4(3m — it) (26) 
ctr ar 

which lead to a linear dependence of $0 and Xq on the slow time variable r. 

As is well known the NLS eq. is invariant under Galilei transformation. The 
same is true also for the 1-NLS eq., so without any lost of generality we can solve 
(PH) in the reference frame where the soliton is at rest (u = 0) and then through 
the Galilei transformation 

X = X' + 2uV, T = T', 8% = 5y> e V™x'+2iu*T>) ( 27 ) 

we can find the solution in the laboratory system. The one-soliton solution in its 
own reference frame is given by 

2iv 2 T 

* = 2w — : -. (28) 

cosh 2vX K ' 

We are looking for solutions of the form 

M = ie ^y{X) (29) 

with Y(X) a real function. With 2vX — > X and introducing the new variable 
p = tanh X the equation satisfied by Y becomes 

i (<* - K ) + < 2 - Th )Y = 48 * V ^ (30) 

In the lhs we recognize the equation for the associated Legendre polynomials. 
The two linear independent solutions are P\ and Q\ 

1, 1+P , P 



pI = -^-p\ & = -^-p\2 ln T^p- + T^)- (31) 

We write the solution of the nonhomogeneous equation fl3"0|) as a linear superpo- 
sition 

Y(p) = a(p)P^p) + (3(p)Q 1 1 (p), (32) 

where a(p), f3{p) are coefficients to be determined. A general way to find them is 
to consider 

Then using the Wronskian expression it is easily found that 

f(p) = -2Av 2 P Jl^. (34) 



Then the equations for a and (3 are easily integrated giving 

a(p) = -3 V 2 (1 - p 2 fln\±^ + 6v 2 p(l + p 2 ), /3(p) = Qv 2 (l - p 2 ). (35) 

1 - p 

In integrating the equation for f3 an integration constant was determined from 
the condition /3(±1) = 0. This ensure us to have Y(p) finite everywhere for 
p G [—1, +1]. Now introducing Q3l)| ) in (^) a very simple form for Y(p) is found, 
namely 

Y(p) = -I2v 2 p^l - ?. (36) 

As Y(p) is zero both in the origin and at p = ±1 (X — > ±oo) it has a maximum 
at p = ±-^, which transformed in the X variable gives X = ^ln(l + \/2). This 
maximum is similar with that found in numerical simulations and other theo- 
retical treatments of the similar ZDP problem of pulse propagation in nonlinear 
optical fibers [0-0. 

4. Conclusions 

Solitonic type excitations in a Davydov model are investigated. The multiple 
scales method is used to study the space-time modulation of the amplitude. As 
expected the dominant amplitude A\ is satisfying a NLS eq. In the next order 
the second amplitude is given by the solution of a nonhomogeneous linearized 
NLS equation. Possible secular behaviours are generated if symmetries of NLS eq. 
are identified in the rhs of this equation. They are eliminated if the NLS solution 
satisfies also the next eq. in the NLS hierarchy (a complex mKdV equation). 
The ZDP case is also investigated using the same method of multiple scales, but 
with another definition of the slow spatial variable. The dominant amplitude 
satisfies a modified NLS equation containing a third order derivative term. As its 
coefficient is a small quantity a perturbational approach is used. The case of the 
one-soliton solution of the NLS eq. is fully solved. Possible secular behavior is 
eliminated if the one-soliton solution satisfies also the next equation in the NLS 
hierarchy. A linear r 3 dependence of the initial phase and position of the soliton 
is determined. The remaining equation, by a suitable transformation, is reduced 
to the equation for associated Legendre polynomials. The complete solution is 
found if the NLS soliton is at rest. By a Galilei transformation the solution for 
moving soliton is obtained. More complicated situations are under investigation. 
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